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The zinc-blende GaAs-based III-V diluted magnetic semiconductors (DMS)are studied in the 
coherent potential approximation(CPA). In this work, we use the exact Hilbert transformation of 
the face-centered cubic(fcc) density of states (DOS), which is different from the usual semi-circle 
density of states employed in our previous work. Our calculated relation of ground-state energy 
and impurity magnetization shows that ferromagnetism is always favorable at low temperatures. 
For very weak Kondo coupling, the density of states (DOS) of the host semiconductor is modified 
slightly. Impurity band can be generated at the host band bottom only when Kondo coupling is 
strong enough. Using Weiss molecular theory, we predict a nonlinear relation of Curie temperature 
with respect to Kondo coupling, as is different from the conclusion of our previous calculations based 
on semicircle DOS. The agreement of our calculated Tc with measured values is convincing. 
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The ferromagnetism of DMS of Ill-V-type is not well 
understood. To explain ferromagnetism in DMS, var- 
ious models and approaches have been proposed0-[^. 
Though the models differ from each other in details, they 
all agree that the coupling between the carriers and local 
spins is of fundamental importance. An issue of debate, 
however, is how the exchange between localized spins is 
induced by the carriers. One model for this induced ex- 
change is the Rudernian-Kitttel-Kasuya-Yosida (RKKY) 
interactional Another version which results in con- 
clusion equivalent to RKKY is the Zener model|4j which 
uses the fact that the valence holes are on p-orbitals. A 
third model is the double-exchange (DE) mechanismQ. 
But this model is inconsistent with the charge-transfer 
propertiesy. Though RKKY can give a Curie temper- 
ature in agreement with experiment, some argue that 
the RKKY model breaks down here 5j 8j because the lo- 
cal coupling between the carrier and the impurity spin is 
much larger than the Fermi energy and can not be treated 
perturbatively. In dealing with the effect of the localized 
spins, a key issue is whether or not randomness should be 
taken into consideration. The above models are all mean 
field approximations (MFA) which assume homogeneity 
and neglect randomness. But DMSs are disordered sys- 
tems with positional disorder of Mn impurities. As con- 
cluded in 2] 9] , disorder has a substantial influence upon 
carrier magnetic susceptibility. Hence, any flrst principle 
consideration should take into account the randomness 
of the impurities. 

A classic method of dealing with randomness is the co- 
herent potential approxiniation(CPA)fl^ which has been 
applied to DMS llj- 13] . Basing on the formalism of[l^ 
and assuming very large local spin S while keeping the 
product IS constant (where / is the Kondo -like inter- 
action), 1^ obtained the density of states and the re- 
lation between Curie temperature and the doping con- 



centration. Using the averaged carrier Green's func- 
tion, |0| arrived at the conclusion that the local cou- 
pling between the carrier and the impurity spins must 
be intermediate in order to acquire ferromagnetism. In 
our previous workfis). we used the formalism of CPA in 
[l^ to study the ground-state properties of III-V DMS. 
The DOS for the undoped crystal we used is the semi- 
circle DOS. In contrast to we kept S = 5/2 and 
treated the impurity spins fully quantum mechanically. 
Though it is mostly accepted that in III-V typ e DMS 
the effective spin of valence holes is we de- 

scribed the holes as spin-1/2 fermions. It is usually be- 
lieved that such a description can still catch the essential 
physics. Because of spin-orbit interaction, the p-orbitals 
are spilt (with split-off ~ 0.34eV) into a spin-3/2 niul- 
tiplet and a spin-1/2 multiplet|2^. Using spherical ap- 
proximation, the kinetic energy of the Luttinger-Kohn 
Hamiltonian for the spin-3/2 multiplet takes the form 
X]^(S^k^/2TOp)C]{j.^Ck^ near the valence top after diago- 
nalization, where = ruh ~ 0.5m for fj, — ±3/2 and 
iTT-p, = iTT-i — 0.07m for fi — ±1/2 (m is the effective 
mass of a free hole ) . The interaction between the spin of 
holes and local 5/2-spins now takes a k-dependcnt form 
Ei,k,k' Si-c|^J(k,k')ck' exp(-i(k-k')-Ri)_19]. Since the 
DOS for parabolic band is g{e) = il/2n^h^){2mf/^^, 
we have the ratio of DOS gh{£) / gi{e) — 19 for heavy 
holes and light holes, i.e., most of valence holes are heavy 
holes. Therefore, it is a valid approximation to con- 
sider only heavy holes. Moreover, since the hole den- 
sity is very small, the Fermi wave vector is supposed to 
be very small and the it is thus a reasonable approxi- 
mation to consider J(k, k') ~ J(0,0) in the interaction 
term. And this leads to the usual assumption that the 
carriers are shallow holes and the coupling of the shallow 
holes to the Mn^+ can be described by local Kondo in- 
teraction between 5/2-spins and 1/2-spins. The system 
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is highly compensated |23|-[25j and the carriers are holes 
originating from randomly distributed Mn. There are dif- 
ferent kinds of randomness, e.g. substitutional random- 
ness, interstitial randomness, antisite randomness, and 
directional randomness of impurity spin. It is commonly 
agreed now that interstitial Mn atoms and antisite As 
only reduce the hole densities and do not affect conduc- 
tion of holes significantly. Therefore we consider only 
two kinds of randomness, i.e., the random substitution 
of the Mn atoms and the random direction of the impu- 
rity spins. 

Since the DOS for undoped crystal is much in- 
volved in the CPA calculations and the real structure of 
Gai-xMuj,, where the doping concentration x varies from 
0.015 to 0.08 in region of interest for ferromagnetismfjfij . 
is of zinc-blende structure, we here use the method devel- 
oped in our previous work to study DMS Gai_a,Mn^As, 
taking into consideration the real fee lattice structure. In 
as-grown samples of , ferromagnetism was realized at 
Curie temperature of 76K at a; = 0.05. The hole density 
p varies from 5.6%a; ~ 93%x. We use the tight-binding 
model Hamiltonian 



and 



H 



(1) 



where u\ depends on whether i is a Ga or Mn site. 
For Ga-site Ui ~ uf = Eq Cj^Qcr , and for Mn-sitc 
Mi = = Em X^cr Aa^'^'^ + "^KSi • Si. Si is the local spin 
of Mn at site i, s = (1/2)cJ.Tctct'Cct' is the spin of a hole 
where cJ.(cct) is the creation(annihilation) operators for 
holes , spin indices a, a' =t, i, and r = (n, T2, ra) are the 
three usual Pauli matrices. Em and Eq are the on-site 
energies for Ga and Mn and are assumed constant. The 
hopping energy iij = i if i,j are nearest neighbors and 
zero otherwise and Jk > is the local Kondo coupling. 
According to the general scheme of CPA, the virtual un- 
perturbed Hamiltonian is 



Vtk + 5]^(e))cLck^ 



(2) 



where e is the Fourier frequency variable, ^^(e) is the 
CPA self-energy to be determined self-consistently and 
tk is the Fourier transformation of tij . Then the relative 
perturbation V is given by 



(3) 



where Vi = vf ^ Y^A^G - ^a)c\^Ci„ for Ga and 

^'i = < = Ea(^M - ^MaCi<y + J^^i ■ Si for Mn. The 
reference Green's function is (icr|^(e)|jcr') = {0\cia{e — 



Jif) Cj^/|0) where |0) is the vacuum state of the c- 
operators, and the associated t- matrices are t9 — Wi'^/(1— 
^v9)^tf^ = t;M/(l - Mvf). So the CPA equation and 
DOS are given by 







(4) 



ga{e) = --ImF^(£) 

TT 



(5) 



where (• • •)spin denotes average over the configurations 
of impurity spins and Fcr{e) is the Hilbert transforma- 
tion of DOS, F^{e) = {ai\^\ia) = (l/iV)Ek[V(£ ^ 
tk — So-)], where N is the number of lattice sites. For 
any /(S'^), the spin average is given by (/(S"=))spin = 



f{S^)/ J2s'=-S ^ where A is determined 



by the condition {S^ 



' spin 



m, m is the given magnetiza- 



tion of the impurity spins. In our single-particle CPA, the 
Callen-Shtrikman relation[2^ that tells there is a one-to- 
one correspondence between m and ((S'^)")spin for n > 1 
applies . These explicit forms of CPA equations Q are 
given by in another context. To solve CPA equa- 
tions, we need to know Jelitto0| and Iwata 
29] first gave a rigorous calculation of the k-integrated 
Green's function for fee tight-binding Hamitonian. The 
expression entails complete elliptic integral and the an- 
alytic properties are very complicated. Morita 30] et al 
made further discussions later. A much improved ex- 
pression of the k-integrated Green's function is obtained 
quite recently 01 . In unit of w/A where w is the band 
width and defining 



1., Ay^VT^ - (1 + e)VT^V-3 - s, 
<^ = ] 

then F(o)(£) can be written as 
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(6) 



7r2(l-e)3/2 



(\/-£-3-2\/^)/C2(k) (7) 



where /C(k) is the complete elliptic integral 



Tr/2 



dO- 



(8) 
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At zero temperature, the carrier density for spin a can 
be expressed as = J^^gai£)de, where ep is the 
Fermi energy, and the total carrier density n = n| -I- n^. 
The total electronic ground-state energy per site is Eg — 

/!^e[5T(e) +5i(e)]'^£- 

In our numerical calculation, we chose w — 4.0 eV 
which is roughly the bandwidth of the heavy holes [s^ 
and set Eq = since we can shift the chemical potential 
without loss of physics. (Note that the effective mass for 
tight-binding spectrum is related to the bandwidth by 
nih = Sli^ / (wa^) where a ~ 5.6532A. This relation gives 
rrih ~ 0.48me for w = 4 eV). FICQshows DOS for some 
given model parameters. The corresponding hole den- 
sity is p ~ 0.442a; for as-grown samples as provided in 
|26j. The left panel shows DOS in the full range for some 
given model parameters; the right panel is an enlarged 
display for the range e = —4 to e = —2. As a check for 
our numerical results, the sum rule ga{z)dz — 1 and 
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FIG. 1: Left panehDOS in the full range for some given model 
parameters; Right panel: enlarged display from e — —4 to 
e — —2, the Fermi energies are for as-grown samples from 

m 



FIG. 3: Theoretical estimation of Curie temperature using 
Weiss molecular theory 
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FIG. 2: Relation of Eg versus m. B: Jk = 0.6; C: Jk = 0.8; 
D: Jk = 1.0; E: Jk = 1.2; F: JK = 1.4; G: Jk = 1.6; H: 
Jk = 1.8. X = 0.05, p = 0.442a;, w/A = 1.0 eV 



the relation + — p are also preserved where n^r is 
calculated from spin-resolved DOS respectively. As can 
be seen that for Jk = 1.2, there is no separate impurity 
band and the DOS is not substantially different from the 
bare DOS in shape. When Kondo coupling Jk = 1.8, 
separate impurity bands show up. Fig.l also shows that 
for TO > 0, there are always more spin-down carriers than 
spin-up carriers, in compliance with the fact that the lo- 
cal p-c? coupling is antiferromagneticj^. 

FIG. 121 shows the relation between ground-state en- 



ergy per site and the impurity magnetization m. In ac- 
tual calculation of ground-state energy, an important is- 
sue is the determination of the Fermi energy which is 
fixed by the integration of an interpolation function of 
the DOS. If the interval 6e is not small enough, the inte- 
gral of the interpolated DOS may vary significantly with 
the choice of Se. To make it stable enough, we choose 
fe = 5 X 10^* to interpolate DOS linearly. Our results 
show that for all the chosen values of model parame- 
ters, the ground-state energy per site always decreases, 
though very slowly, with increase of impurity magnetiza- 
tion. Therefore CPA predicts that at very low temper- 
atures, ferromagnetism is always energetically favorable 
for all the model parameters considered. 

In FIG. 13 we show the dependence of Curie temper- 
ature on the model parameters and the doping concen- 
tration. As proposed in our previous work [l5j, Weiss 
molecular field theory can be employed to calculate the 
Curie temperature. Given m, one can calculate DOS and 
then (sz). So one can establish relation (s^) = {sz){m). 
On the other hand , given (s^), each impurity spin feels 
an effective field Jk(sz) and thus we have to = SBs{(3h) 
with h = — Jk{s z) ("m) , f3 = l/k^T and Bs{x) the con- 
ventional Brillouin function. For very small m, we have 
{sz) — —Am with A > and we have f3h ~ (iJ^iAm . So 
Bs{(3h) ~ (5 -I- l)/3/i/3 and thus the Curie temperature 
can be estimated by k^Tc — JkS{S + 1)^/3 . For small 
TO, letting (z) — F{z) + 'ip{z)m, (z) — F{z) — ip{z)m 
where F(z) is the paramagnetic solution, we have 



A{Jk,Em,x,(3) 



1 



de lmtp{e) 



(9) 



where we have ignored the /3-dependence. As our nu- 
merical results (not shown here) indicate, the chemical 
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potential is very close to the zero temperature Fermi en- 
ergy in a wide range of temperatures (/3 > 50), showing 
that the Fermi function can be approximated by the zero 
temperature step function. The left panel in FIG. 3 shows 
the relation of Tq versus Jk for different values of Em- 
Unlike the result from semicircle DOS the curves 
exhibit a nonlinear relation. The right panel shows the 
dependence of Curie temperature on the doping concen- 
tration for Jk = 2.03, Em ^ and Jk = 1.2, Em = -3.25 
for the as-grown samples in [2^. These two sets of pa- 
rameter values are chosen so that at a; = 0.05, Tc = 76 
K. It is seen that our CPA -|- Weiss approach for Tq does 
reveal the observed behavior. Our result show that both 
two sets of values Jk,Em can explain the measured re- 
sults. 

To conclude, we summarize our results here. Using 
CPA and treating the impurity spins fully quantum me- 
chanically, we have calculated the ground-state energies 
of GaAs-based III-V DMS for a wide range of model pa- 
rameters. The results show that ferromagnetism is al- 
ways preferable at low temperatures . With the help of 
Weiss molecular theory of ferromagnetism, we obtained 
a nonlinear relation of Curie temperature with respect 
to Kondo coupling. Our theoretical prediction of rela- 
tion Tc ~ a; does fit quite well the measured results for 
as-grown samples |^ . The deviations of our calculated 
results from the experimental values is due to the sensi- 
tivity of Tc to the hole density or the ratio p/x, which 
is difficult to be measured precisely. Though the value 
of Jk = 1.2 eV is more generally accepted than other 
values, the actual value of Jk is still an issue of debate. 
Our calculation allows adjustment of model parameters 
Jk, Em and the main features of numerical results for dif- 
ferent model parameter values remain similar. Finally, a 
number of complicated models with complicated approx- 
imations have been proposed to calculate Curie temper- 
atures (e.g.jl^)- Yet, we believe that CPA, a method 
having been tested for decades, is still an effective way 
to understand ferromagnetism in DMS. 
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